Entanglement in a model for Hawking radiation: An Application of Quadratic Algebras 
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ABSTRACT 

Quadratic polynomially deformed su{l, 1) and su(2) algebras are utilised in model Hamiltonians to show 
how the gravitational system consisting of a black hole, infalling radiation and outgoing (Hawking) radiation 
can be solved exactly. The models allow us to study the long-time behaviour of the black hole and its out- 
going modes. In particular, we calculate the bipartite entanglement entropies of subsystems consisting of a) 
infalling plus outgoing modes and b) black hole modes plus the infalling modes.using the Janus-faced nature 
of the model.The long-time behaviour also gives us glimpses of modifications in the character of Hawking 
radiation. Lastly, we study the phenomenon of superradiance in our model in analogy with atomic Dicke 
superradiance. 
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Classically, black-holes, being the final state in the evolution of massive stars, cloak their singularities 
lO from the rest of the spacetime through event horizons and have surprisingly simple properties. Ob- 
servers who live in the region exterior to the event horizon cannot observe the singularity inside the 
event horizon while those that have travelled inside cannot communicate with observers outside the 
qh event horizon. Black hole solutions, which are "eternal", can be described at most by three parame- 
ters, their mass, angular momentum and electric charge. A rotating black hole can lose energy through 
W) an energy extraction process (Penrose process- which from a thermodynamic perspective has been 
called incoming wave amplification, or superradiance with associated "grey body" factors for the black 
hole HI), classically, black holes cannot decrease but can only increase the area of their event horizons 
^ (through accretion of matter) . 

A dramatic shift in the black-hole paradigm occured when Hawking, using earlier results that a clas- 
sical black hole cannot decrease the area of its event horizon, made a remarkable prediction that from a 
quantum mechanical prespective, a black hole can emit radiation which an observer at infinity sees as 
being thermal. It was suggested that the process itself is through quantum fluctuations in the intense 
gravitational fields near the event horizon of a black hole. Quantum fluctuations create pairs of parti- 
cles, which are split apart, one falling into the black hole and the other radiating away in such a manner 
that the black hole appears to be emitting thermal radiation at late times |2|. The energy required for 
^ this radiation is believed to come from the mass of the black hole, which eventually disappears. Hence, 
from a quantum mechanical viewpoint, the area of the event horizon decreases (with a corresponding 
decrease in its mass) until the black hole eventually ceases to exist when its mass goes to zero. These 
studies of quantum fields in interaction with black holes has led to extraordinary links between black 
holes, thermodynamics and information theory |3 |. Because of this interaction with quantum fields, a 
black hole can be associated with temperature related to the surface gravity of the horizon. In its sim- 
plest form for uncharged, non-rotating black holes, T oc 1/M, where M is the mass of the black hole. 
This inverse relationship between the mass of the black hole and its temperature gives us an impor- 
tant clue to its evolution. In particular, Hawking has shown that the black hole behaves as a black body 
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at this temperature and correspondingly, an observer at infinity sees thermal radiation from the black 
hole. Energy conservation would then imply that as the black hole radiates energy, its mass decreases 
with a consequent increase in temperature. It is clearly important to understand this later stage evolu- 
tion of the black hole as it eventually disappears. More recently it has been suggested that if a quantum 
theory of gravity is to explain Hawking radiation, then there would be a clear violation of quantum me- 
chanics |4|. This is caused by the evolution of pure states (states falling into black hole) into mixed 
states (thermal radiation from a black hole). Notwithstanding a string field theory explanation (5j[6], 
it is important to analyse the complete evolution of a black hole through to its disappearance, while 
examining the states radiated out to infinity. This is only possible through simple models, where, back 
reaction effects can be included in later stages of the evolution where they would dominate. Of course, 
we would like to answer the question of the information paradox being resolved through an audit of the 
complete radiation released out to infinity by the black hole through its lifetime. Since it is very difficult 
to actually test such cosmological predictions, intra-disciplinary links in physical systems are especially 
important and allow analogue laboratory models to serve as important tools in their understanding. 

Laboratory analogues for simulating Hawking radiation are now blossoming, with the advent of new 
techniques in condensed matter physics and quantum optics 1 7 1 . Such models, of course, are at the most 
toy models, and examples are acoustic models in moving fluids ( in which supersonic fluid flow gener- 
ate an acoustic analogue of a "black hole" (8], and the presence of phono nic Hawking radiation from 
the acoustic horizon is derived), optical models which use slow light passing through a Bose Einstein 
condensate that is itself spinning faster than the local speed of light within to create a vortex capable 
of trapping the light behind an event horizon to produce an optical black hole[9|. Recently using pico 
second pulsed lasers passing through glass, the conditions for production of a disturbance zone in an 
optical black hole have been produced in the laboratory. This disturbance zone exhibits the emission 
of analogue hawking radiation for an "optical black hole" 1 10 1. There are many other "analogue mod- 
els" that may shed new light on perplexing theoretical questions concerning the information paradox 
problem. The information flow is in principle bi-directional and sometimes insights developed within 
the context of general relativity can be used to understand aspects of the analogue model. We seek to 
analyse a model that keeps the salient features with respect to just the information loss aspect, enabling 
us to see whether the thermal nature of the radiation is modified when the black hole has evolved to 
a Planckian size and is therefore amenable to quantization. We consider a first quantized model as an 
illustrative example in order to examine if there is any non-thermal contribution to the Hawking type 
radiation arising from the quantization of the analogue "black hole". 

We consider two such models. The first is the model introduced by Nation and Blencowe 1 1 1 1 . The sec- 
ond is the Tavis-Cummings version of the Dicke model|[T2]. The two mode parametric amplifier with 
a bilinear Hamiltonian has traditionally been used as a toy model for hawking radiation, since it pro- 
duces two mode squeezed states, which can be mapped on to the incoming and outgoing states across 
the black hole horizon and thus simulate Hawking radiation. However, if we want to consider the effect 
of back reaction, we have to include the coupling of the quantized "black hole" mode. This is achieved 
by using a trilinear Hamiltonian [13] . Traditionally in quantum optics such a Hamiltonian is used to 
describe a quantum parametric amplifier and frequency converter with a quantized pump field. The 
outgoing modes are identified with the signal mode, the incoming modes are identified with the idler 
modes and the gravitational field of the black hole is represented by the pump HU . It can be shown that 
if the pump mode is treated classically, the signal mode (corresponds to the outgoing particle modes) 
and the idler mode( corresponding to ingoing particle modes) form a two-mode squeezed state, and if 
we trace over the idler modes, the outgoing radiation can be mapped onto Hawking radiation. In the 
trilinear version, the pump mode is also quantized . This, thus, corresponds to the quantization of the 
black hole. A nice pictorial analogue can be found in HT][T5]. For our purposes in this paper, our no- 
tations are such that mode b is used to represent outgoing particles (Hawking radiation), c is used to 
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represent the infalling particles and a is used to represent the quantized black hole modes. Essentially, 
the modes b and c represent the pair creation process at the event horizon of the black hole. 

Nation and Blencowe have studied the entanglement properties of this system using a short time 
approximation! 11 1. Our interest is an analysis of the complete history of the radiation and the short 
time approximation is insufficient for our purpose. We are able to do a comprehensive study by ex- 
hibiting a Quadratic Polynomial Algebra (a polynomial deformation of a Lie algebra) of the trilinear 
Hamiltonian in two different ways. One of which is similar to that of HU, while the other maps the 
system onto an interacting atom-radiation system called the Tavis-Cummings model. This Janus-like 
nature of the mapping, enables us to study many different aspects of Hawking radiation in analogue 
black-hole systems. 

n. CLASSICAL BLACK HOLE 

The general Hamiltonian in the model system of interaction between the particle and black hole 
modes is given in the trilinear form 

J€ = co a a^a + co b b^b + co c c^c + k ab^c* + k* a} be . (1) 

For ease of calculation and comparison with the results of [TT), we define k = ih% to get 



J€ ' -co a a* a + co b b^b + co c c 1 c+ ih%ab^ V - ih%* a* be. (2) 

Here a, b, and c are the boson annihilation operators of the quantized black hole modes, the outgoing 
modes and the incoming radiation modes respectively. Energy conservation requires co a = co b + co c . In 
order to show the entanglement between the modes b and c, we first treat the black hole mode a as 
classical (c-number A) to give the Hamiltonian 

J€ = a) b b^b + coc c + c + ih%Ab ] c f - ih%* Abe . (3) 

This Hamiltonian can be written in terms of the generators of an su{l, 1) algebra 

K + = b f c\ K- = bc, K Q = - (4) 

with Lo b = co c = co, as 

je = lo{2Kq - 1) + ih%AK+ - ihx* AK- . (5) 

The corresponding Casimir is 

<X = K%-^— 2 + - = . (6) 

The representations are labelled by the eigenvalues of the Casimir operator and the number operator 
for the outgoing modes/^: 

Jf\k,nb) = k(k-l)\k,nb) (7) 
The action of the su(l, 1) generators on the basis | k, n b ) is 

K \k, n b ) = [k+ rib) | k, n b ) , (8a) 

K+ \k, n b ) = \/{n b + \){nb + 2k) \k, n b + l) , (8b) 

K- \k, n b ) = \J n b {n b + 2k- 1) \k,n b - l) (8c) 
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A convenient basis to use is the two-mode oscillator basis 



\n b ,n c )= - — - |0,0> . (9) 



n b \ n c \ 

We set n c - n b = q, with q being called the degeneracy parameter and takes values q = 0, 1, • • • . It is easy 
to go from one basis to another by noting k = |(1 + q). Since the Hamiltonian is linear in the generators 
of the Lie algebra of the SU(l, 1) group, the unitary evolution due to the interaction is given by 
the coherent state relevant to the discrete series representation of SU(l, 1). 

\y{t)) = exp[P(t)K + -p*(t)K-}\k,0). (10) 

Where /3(f) = xAt and the kinetic part contributes an overall phase factor. Using the disentanglement 
theorem (Zassenhaus formula) for su(l,l) we get 

exp[p(t)K + -p*(t)K-] =exp[ r 'K-]exp[i 1 Ko]exp[rK + ], (11) 

where, y = tanh(|/3(f)|)e^, rj - 2lncosh(|/3|), and y' = -y* . This allows us to write the time evolved wave 
function as: 



. °° Y{2k+n b ) , 



In the two mode oscillator basis this becomes: 



rib 

For the special case q = (overall neutrality), 



TO>> = (i-rr +11 '%/^r"»K«>+<7> 

o v n b- q- 



(13) 



|^U)) = ^i-r 2 £ r" 6 1 (i4) 

The density matrix is 

Cm b +q,m b \ n b> ™b + Wfc + g, m fo | , (15) 

nb> m b 

where 



The entanglement properties are calculated by taking the partial transpose of p : 

p(f) = £ C nb , mb+? C* b+ ^ m Jni,,ra fc + <7)(ra & + 4,rai,| . (17) 

/If,, TO;, 

When n b = n c = n ,{q = 0), 

Ppr(f) = £ c «,™ C to,« \n,m)(m,n\ (18) 

h,to 

Where C w>m = tanh m+n (^ylt)sech(jylt). The Eigenvalues can be found by looking at the terms of the 
matrix fT6l 

c n, n C* m m \n,m)(m,n\ + C mim C* n \m,n)(n,m\, (19) 
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FIG. 1: P T i b vs. rif, for the outgoing modes 

and are 

A nn = tanh 2 "(jA?)sech 2 (xAf); n = m (20) 

X nm = ±t<mh m+n (xAt)sech 2 (xAt); n^m (21) 

Since the eigenvalues can be negative , the state is entangled, by the Peres-Horodecki citerion. The 
probability distribution of the outgoing modes is given by P nb = seek 2 {Ax) tanh 2 " b (Ar) and a plot of 
P nb vs. rib for different values of the scaled time t is shown fig. |[TJ. 



The Entropy of entanglement is given by 

S = -Tr„plogp = cosh 2 (j At) log cosh 2 {x At ) "log sinh 2 (^ylf) (22) 

Figure |2j shows how the entanglement entropy increases with the scaled time t and it is clear that as 
t — ► oo the states b and c are maximally entangled. 



To show the thermal nature of the outgoing radiation, we define the temperature T = 
[H(i) c l2kB ln(coth(Ar)] allows us to write the entropy in a more familiar form: 

Jan. hco Jbl , 

S = -log(l-e^)-— (1-e^r 1 (23) 

K B 1 

The probabilty distribution of the outgoing modes, P„ b (l/r) = e~ 2IT [e llT - 1] is plotted in fig. |3| as a 
function of 1 / T, showing that the outgoing modes have a thermal distribution in the form predicted by 
Hawking for black holes. 



We therefore see that in the long time limit there is maximal entanglement between the incoming and 
outgoing modes. We also see that a thermal distribution results for the outgoing radiation, as expected, 
when the black hole is considered to be a classical object. In the next section we quantize the black hole 
modes and reexamine the thermal and entanglement properties. 
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FIG. 2: Entropy of the outgoing modes corresponding to the case of a classical black-hole. 
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FIG. 3: Planckian distribution of the outgoing radiation 



III. QUANTIZED BLACK-HOLE 

We now come to the advantage that analogue models provide. In the absence of a quantum theory 
of gravity, it is only through models such as the one we are studying, that we can actually examine the 
result of quantizing the "black hole" modes. Such effects have only been studied in the short time limit 
fTTl . which clearly is not sufficient to determine effects at infinity as required by black hole physics. 
We would like to see if the thermal properties and the entanglement properties in the long time limit 
change when the black hole is quantized. In this section, we shall show how quadratic algebras can be 
used to solve this particular problem. 

The full Hamiltonian with three quantized modes can be written as in eqn.Q, 

J€ = (i) a (a i a + K ) + ih%aK+ - ih%* a^K- , (24) 
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by choosing w/, = a> c = co a /2. 
Denning 



the Hamiltonian Eq. 124 1 takes the form 

je = 2o) a &+ihxQ+-ihx*Q-. (26) 

These operators Q+, Q_ and Q satisfy a quadratic algebra, which is a special case of Polynomially de- 
formed algebras, introduced first by Sklyanin [17], in the context of braid groups and modified to be 
used in quantum optics by Karassiov and Klimov (T8) . This quadratic algebra is given by 

[Qo,Q±] = ±Q± 

[Q + , Q_] = 3 + {2& - 1) Qo - {X + 5£{S£ + 1)), (27) 

where, 5£ is the central element of the quadratic algebra and JT, defined in equation!^, is the Casimir 
for the su{l, 1) algebra. The Casimir Operator for this algebra is 

c € = Q+Q- + Ql + {££-2) Ql-{Jf + S£ 2 + 2S£-l) Q + [<X + ££{££ + (28) 

This is a special form of the generic su{\, 1) quadratically deformed polynomial algebra 

[Q+, Q_] = a Ql + p Qo + ( = g(Qo) - g(Qo - D = /(Qo) (29) 

where the function, g is called the structure function and is defined upto the addition of a constant as, 

a o (a + B) ? (a + 36 + 60 

g(Qo) = -Qo 3 + — 7-^Qo + ^ — -Qo ■ (30) 

o 2 b 

The Casimir expressed in terms of the structure function is 

^ = Q_Q + + g(Qo) = Q + Q- + g(Q -l). (31) 



The quadratic algebra (|27} corresponds to the choice a = 3,/3 = {££ - 1), and £ = - \J£ + ££{££ +\)\. The 
quadratic algebra belongs to a class of algebras similar to the enveloping algebras of 5/(2, C) and has a 
notion of highest weight modules and every finite dimensional module is semi-simple HU . 

Thus, the basis | k, £, rib) of a unitary representation of the quadratic algebra is finite dimensional and 
characterized by the quantum numbers k and £ [20] . The action of the generators is given by 

££ \k,e,n b ) = e\kj,n b ) , (32a) 
Qo I k, £, n b ) = (k + n b - £ ) | k, £, n b ) , (32b) 

Q+ \k,e,n b ) = \/{n b + l){n b + 2k)(2£-k-n b ) \k,£,n b + l) , (32c) 

Q-\kJ,n b ) = \/n b (n b + 2k-l){2£-k-n b + l)\k,£,n b -l) . (32d) 

With n b = 0,1,2,- ■■ ,(2£-k). 

For our purposes it is better to redefine the states in terms of the following constants, n c - n b = q = 
2k-l , p = n c + n a and p - q = 2£ - k. With the new states, the action of the generators can be written 
down as follows: 

i , (4n b + 3q-2p + l\ . . 
Qo\p- q- n b ,n b + q,n b ) = ^ j \p - q - n b ,n b + q,n b ) , (33a) 

Q+ \p-q- n b ,n b + q,n b ) = \J{p -q-n b )(n b + q + l)(n b + l)\p-q-n b -l,n b + q + l,n b + l) , (33b) 
Q-\p - q - n b ,n b + q, n b ) = \J{p - q - n b + l)(n b + q)n b \p- q - n b + 1, n b + q- 1, n b - l) . (33c) 



8 



From the structure of the Hamiltonian, it is evident that the incoming and outgoing particle modes 
are being created through the annihilation of the black hole modes. Therefore we have the following 
conditions imposed by particle number conservation: n c = n b , n a = Wc 2 " b , hence, q = and p = 2n b . 
We start with a state in which all the energy modes are concentrated in the black hole. Thus, 

i . (4n b -2p + l\ | . 
Qo\p- n b ,n b ,n b ) = \\p-n b ,n b ,n b ) , (34a) 



Q+ \p- n b ,n b ,n b ) = \J{p - n b ){n b + \){n b + \) \p- n b - l,n b + l,n b + l) , (34b) 
Q_ \ p-n b ,n b ,n b ) = \J{p-n b + \)n b n b \p - n b + l,n b - l,n b - l) . (34c) 

With the inclusion of the quantized black hole modes, the Hamiltonian can be written as a linear 
combination of the generators of the quadratic algebra. Since S£ is the central element of the quadratic 
algebra , the time evolved state of the Hamiltonian is the Perelemov coherent state (PCS) of the quadratic 
algebra: 

| Y(T)) = expT(Q + -Q_) |p,0,0) (35) 

where t is the dimensionless time % t. The ground state deserves some comments. The initial state is 
not completely devoid of energy, this is indicated by the degeneracy factor p ^ 0, in other words, in the 
initial state all the energy is concentrated in the BH. 

As it stands the exponential can not be disentangled due to the nonlinear terms occurring in the 
commutation relation. We have used deformed algebra techniques to develop a Zassenhaus formula 
for quadratic algebras to calculate the coherent states in (T4) . For the sake of completeness, we will 
discuss this construction in some detail. 

The first step is to map the quadratic algebra to su{l, 1) or su{2) as the case may be, this process is 
called the linearization of the algebra. This technique relies on the crucial step of constructing a new 
annihilation operator Q_ such that it satisfies the following commutation relations [Q+,Q_] = -2AQ 
and [Qo, Q-] = -2Q_. When A = ±1 the algebra is either su(l, 1) or su (2) respectively. To find this new 
operator Q_ we make an ansatz 

Q-=F{ c e,Qo)Q-, (36) 

where F is chosen to be a function of the Casimir and the diagonal generator. This function can be 
solved for by substituting this ansatz in the linearized commutator, i.e., 

Q+ F(^,Q ) Q--F(&,Q ) Q- Q+ = -2AQ 
F(^,Qo-l)Q + Q--F(^,Qo)Q_Q + = -2AQ . (37) 



In the present case it is 



Q {Q + l)X + £ 

F{^, Qo) = • (38) 

^-g(Qo) 



Where e is an arbitrary constant that can be fixed by noting that the ground state of Q_ and Q_ is same. 
Note that in the present study A = 1. Details of this construction can be found in fT4l . 
Now that we have a new lowering operator, the PCS is easy to construct and is given by 

|¥(T)) = ex P T(Q + -Q_)|p,0,0). (39) 



Note that Eq. (35 is equivalent to Eq. ( |39| , since in both cases the annihilation operator acting on 



the ground state gives zero and the diagonal operator gives trivial phase and can be absorbed into the 
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normalization constant. Thus the operator responsible for the time evolution of the state is Q+. Having 



said that, one must bear in mind that the evolution given by Eq. 139 is exact and valid for even large 
times unlike the evolution advocated by IfTTI . 



Since the algebra is linearized we can apply the disentanglement formula 1 1 1 1 to Eq. J39) to get 



Kt(t)) = — — exp(yQ+)|p,0,0), (40) 
VN(t) 

where N(r) is a normalization constant and y = tanh(r); r = %t- Using Eq. f34b| l, the CS is 



\w : yM)= - 1 V / Hp + D r nb \p-n b ,n b ,n c ) , (41) 

and the normalization constant is 

N(j) = exp[coth(r) 2 ] tanh 2p (r) T{p + l,coth(r) 2 ) . (42) 

T(a, b) is the reduced gamma function. For the sake of completeness the above CS in the | k, £, n b ) basis 
is 



, > 1 2£ ~ k T{2k+n b )T(2£-k + l) _ . „ . 

Y,r(T)) = -== Y \ „ = r b \k,e,n h ). (43) 

1 ' 1 \/r(2A;)r(n fo + l)r(2/-A:-n fo + l) ^ 

The normalization constant for the above CS is given by the Tricomi confluent hypergeometric function 
U{a, b, x), 

N{r) =j~ Ak U[2k,2£ + k +1,7" 2 ) . (44) 

The state of the system consisting of black hole (a)+incoming modes (c) + outgo ingmodes (b) is a 
pure state and therefore the total entropy is zero. Hence the marginal entropies of the black hole and 
the particle-antiparticle subsystems are equal: S a = S bc . From the Araki-Lieb |21| theorem, \S b - S c \ > 
Sbc = S a >\S b + S c \.A quantitative measure of the entanglement between two subsystems is the index of 
correlation I x - y = S x + S y - S xy . In our case the index of correlation J fl _i, c between the BH and Particle- 
antiparticle subsystems is equal to twice the marginal entropy of Black Hole:/ fl _fo c = 2S a . 

The quantum nature of the dynamics induced by the trilinear boson Hamiltonian leads to strong 
entanglement between the black hole and the particle-antiparticle subsystems. This can be seen by 
calculating the reduced density operators of the subsystems. 

The time dependent density operator corresponding to the system Eq. pTl is, 



1 A T{p + l)r nb+mb (r)\p-n b ,n b ,n b )(m b ,m b ,p-m b \ 

PabcM = —— JL , -■ (45) 

Nir) „X ^T{p-n b + \)T{p-m b + \) 

In the present case we have much richer structure as compared to the normal su{\,\) case given 
in Eq.|[T5j. The reduced density operator of the black-hole , obtained by tracing over the particle- 
antiparticle modes, is , 

p a = — f r(P + 1) T* n *\p-n b )lp-n b \. (46) 

Using the reduced density operator p a , the particle number distribution (PND) for the BH modes is 

P(n a ) = (n a \p a \n a ) = — TiP + 1) y 1 ^ 2 ^ . (47) 
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FIG. 4: The particle number distribution for the quantized black hole modes for different scaled times t = %t. 
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FIG. 5: Entropy for the black hole modes with scaled time t = %t- 



which is shown in figure Q. We see clearly that the long time behaviour is markedly different from 
the short time behaviour showing that studies utilising short time approximations would not see the 
behaviour exhibited here. The quadratic polynomial algebra methods we have introduced here allows 
us to see the full range of temporal behaviour. 



The entropy of the black hole is given by S a = - P{n a ) log[P(n fl )] and shown in figure {5}. 



The positive nature of the entropy shows that the black hole modes are always entangled. In the evo- 
lution of the black hole modes, there is a rapid increase in its entropy as it interacts with the particle- 
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FIG. 6: Particle number distribution for the outgoing modes. 



antiparticle modes and it then settles down in the long time limit to an entangled state with lower en- 
tropy, having radiated away its available degrees of freedom. Further evidence for this is obtained below 
where we show the entropy of the outgoing modes. 

The effective modification to the Hawking radiation with the incorporation of quantum entanglement 
with the black hole modes can be found by studying various statistical properties of the outgoing modes 
(rib). The reduced density matrix corresponding to the outgoing modes obtained by summing over the 
black hole and incoming modes, is, 



1 & Tjp + 1) 
Pb N(r) n % r{p-n b + l) 



The PND for the outgoing modes is given by 

P{n b ) = (n b \p b \n b ) 

and is plotted in Fig. |6) 



j 2n » \n b )(n b \ . 



(48) 



N(t) Y{p-n b + \) 



(49) 



The plots for the mean particle numbers for the black hole modes and the outgoing modes are depicted 
in figures and l[8j. 



Finally the entropy for the outgoing modes is S b = -Y.n b P( n b) \og[P(n b )], and is shown in Fig. j9j. As 
in the case of the black hole entropy, we see again that the entropy of the outgoing modes also shows a 
limiting, asymptotic value of a finite non-zero entropy and suggests that there may be a final stable state 
of the black hole still entangled with the outgoing radiation detected at infinity. This clearly warrants an 
in depth study of the gravitational black hole with as much quantization that our existing techniques 
allow. 



Let us compare the graphs of Figs. Q and |7) i.e. PND and mean particle number graphs corresponding 
to the BH modes with those of the outgoing modes shown in Figs. (6) and (8) respectively. It is clear from 
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FIG. 7: Mean particle number of the black hole modes as function of time. 
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FIG. 8: Mean particle number for the outgoing modes as a function of time. 
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FIG. 9: Entropy for the outgoing modes as a function of time. 
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FIG. 10: Number Distribution for the outgoing modes as a function of temperature with quantized BH modes (solid line) for p=6 modes, 
and classical BH (dashed line) 



the graphs that at the initial time all the energy is concentrated in the BH modes. As time increases 
particle-hole pairs are created. Those pairs that are created at the event horizon have very important 
role to play. In the pair that is created one particle goes into the BH (ingoing modes) and the other 
escapes the event horizon and can be detected, with a particle number distribution given by Fig. ([6| . 
Let us re-emphasize that these plots are physically relevant for long time only due to fact that we have 
identified a quadratic algebra structure of the model. In similar fashion to the classical black hole case, 
we can introduce a temperature into the system , and compare the resulting distribution of the outgoing 
modes with the Planck distribution shown earlier in Fig. (j3j . 



From the two figures ( 10 and (jTTj) , we that there is a systematic shift from thermal behaviour of the 
outgoing radiation as the number of quantized black hole modes are increased. This clearly illustrates 
the modification of the Planckian nature of Hawking radiation by the quantization of the black hole 
modes and shows that the outgoing radiation is "superradiant". To sum up, in this section we have 
considered the quantum optical analogue of hawking radiation from a quantized black hole as a para- 
metric amplifier and used the polynomial deformation of a su(l,l) algebra. However, the structure of 
the Hamiltonian is such that it can be viewed in a completely different manner, such that the underlying 
symmetry is a polynomially deformed su(2) algebra, which explains the physics of the "superradiance" 
in the system as will be shown in the next section. 
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Particle Distribution vs Planck Distribution 




FIG. 11: Number Distribution for the outgoing modes as a function of temperature with quantized BH modes (solid line) for p=10 modes, 
and classical BH (dashed line) 



IV. SUPER RADIANCE IN BLACK HOLES AND THE DICKE MODEL ANALOGY 



In black hole physics, superradiance is the phenomenon of wave amplification through scattering off 
a rotating black hole and takes place when the wave frequency satisfies the superradiance condition: 

, J is the angular momentum, M ,the mass and r#, the horizon 



JIM 



ir 2 H +(JIM) 2 ) 



< co < mco+ where to+ 
radius of the black hole. 

Superradiance in quantum optics is the cooperative spontaneous emission of photons from a collec- 
tion of atoms. The concept of superradiance can be understood by picturing each atom as a tiny an- 
tenna emitting EM waves. Thermally excited atoms emit light with random intensity oc N. Coherently 
excited atoms have EM field oc N and intensity oc N 2 and leads to Dicke "superradiance" IZZ1 . In this 
section , we reframe the trilinear Hamiltonian model such that the black hole modes and the incoming 
modes combine to form "superradiant bound states" which have an effective spin /. This allows us to 
populate the bound state (atomic) levels in a similar way to that of a rotating black hole with angular 
momentum. We shall see that through the mapping of the Hamiltonian to that of the Dicke model, the 
system allows for superradiance of the Dicke type. By considering a rotating Kerr Black hole as an ex- 
cited state of the static Schwarzchild black hole, Wald and Panangaden |23 | have already observed that 
superradiance is analogous to "Dicke superradiance " in atomic physics . 

Writing the trilinear Hamiltonian ^ in a form analogous to the Dicke model (DM) in the rotating 
wave approximation! 15 1, with the incoming modes and black hole modes forming an atomic system 
(atomic modes) and the outgoing modes interacting with the atomic system, we get: 

J€ DM = (b f b + Jo) + ihxtfj- - ihx*bj+ . (50) 

Jo, J+ form an angular momentum algebra describing the atomic system formed by the black hole and 
incoming modes, which in the Schwinger representation can be written as 

{a* a-rfc) 



h 



cfic , 



J 



ac 



Jo 



(51) 
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Thus, in this case, instead of the b-c system in the earlier section, it is the a - c system that acts like a 
single entity having angular momentum. 
Defining the operators 



P = ^{b i b-J ), P + = tfj-, P- = bJ + , se = hrfb + J ). 



(52) 

z 

the Hamiltonian (|50j becomes 

je DM = 2(o££ + ihxP+-iKx* p - • (53) 

Now the algebra satisfied by P's is 

[P ,P±]=±P±, 

[P + ,P-] = 3P% - (2££ - 1)P + Se{Se + 1)] . (54) 

Comparing with the earlier algebra in terms of the Q's (27) , the linear term in Po is negative, and there- 
fore the algebra is a polynomially deformed su{2) algebra. The basis states of the algebra are product 
states of the su(2) algebra and the oscillator states: 

\j,£,m)-\j,m)\n b ). (55) 

where \n b ) and m)} are the outgoing and su{2) states respectively and rib = 21- m - We construct a 
representation such that £ and !£ are diagonal in the new combined basis 

j \i,e,m) = j{j + i) \j,e,m), se\j,e,m) = e\j,e,m). (56) 

This can be divided into two different cases depending on the value (2£ - j) takes. They are (2£ - j) > 0, 
(2£ - j) < 0. In both these scenarios the representation is that given below 



P \j,£,m) = (£-m)\j,£,m), (57a) 



P + \j,£,m) = J(j + m)U-m+l){2£-m + l)\j,£,m+l) , (57b) 



P_ \j,£,m) = -m){j + m+ \){2£ - m) \j,£, m- 1) . (57c) 

We will now study the time evolution of the state of this model. Proceeding as already shown in the 
case of trilinear system we will construct a new operator P_ such that [P+, P-] = 2P . But for the initial 
state we find it convenient to use the representation wherein instead of m we use n\, making use of the 
central element. Therefore the representation for the algebra now acquires the form 

Po\jJ,n b ) = {n b -£) \j,£,n b ) , (58a) 



P+\j,£,n b ) = JU+2£-n b )(j-2£ + n b + lKn b + l)\j,£,n b + l) , (58b) 



P_ | j, £, n b ) = ^ {j - 2£ + n b ) (j + 2£ - n b + 1) {n b ) \ j, £, n b - 1) . (58c) 
Analogous to the previous case the time evolution of the state is a Perelomov Coherent state for su{2) 

|O(T)) = expT(P + -P-)|j,/,0) > (59) 
where r = %t. Using disentanglement formula for su(2) we get 

|0(t)> = — = exp(rP + )|;,^,0), (60) 
\N (t) 
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where, 7 = tan(r). The resultant time evolved state is : 

ix 1 £ / T(2j -L+l + n b )T(L+l) nu\ „ x 

$(t) =^El — — —Y b \j,t,n b ). (61) 

where L = 2£ + j , N(z) = \ j\ 2L T / ( - L; -2 j; I/I7I 2 ) and T 7 is the Tricomi confluent hypergeometric func- 
tion. The time evolved state can also be written in terms of the atomic modes as: 

fj-L L-j I TU-m + l)T{L+l) w I , 

\<&(T))n M =—r = E A / — f m \j,£,m). (62) 

where f = 7 _1 and ./V flf (t) = pj^ p^y-l+i) ifi ( ~ L; -2 7; |f | 2 ). The density matrix for the state Eq. 161 in 
terms of the field (outgoing) modes is given by 

Y l L 

Pdm = |0(t))<0(t) I = — — E E C(n' b ,L)C*(nl,L)\j,2£-n' b ;n' b )(nl;j,2£-nl\. (63) 

iV( - T - ) »' =0n''=0 

b b 

For the sake of convenience we have abbreviated the summands of Eq. 161 by C{n' b , L) . This can also 
be written in terms of the atomic modes as 

A2{j-L) L-j L-j 

Pdm = I $ (t) )<<D(t) I = ^-— E E C{m',L)C*(m",L)\j,m';2£-m')(2£-m";j,m"\ (64) 

-WflfW fn'——j iYi"——j 

The factors inside the summation in Eq. 162 have been abbreviated by C{m',L). Performing a partial 
trace over the field modes we find the number distribution of the atomic modes as : 

1 , \f\ 2 U-L) TU-m + l)T(L+l) 

Pat = {j,m\pat\j,m) = ^ rfT . -V — \f\ 2m (65) 

N a t(j) T(L-j - m + 1)T(2; -L+ l)T(j + m + 1) 

This is plotted in Fig.fi~2~l 



The entropy for the atomic modes is S at = - Y. m PaAog[P at \, and is plotted in Fig.l 13 



From the plots we see that the distribution of the atomic modes becomes more symmetric with time 
and as the number of modes increases, the entropy increases until it reaches a maximum value corre- 
sponding to a population of all the states and the system is maximally entangled. The particle number 
distribution for the field modes is 

. , . . 1 r(2/-L+l + m,)r(Z,+ l) n„, 

P nb = (n b \p f \n b ) = — — - ; / T "; T —rz \j\ 2nb (66) 

1 1 1 Nf(r) T{n b + \)T{2] - L+\)T{L- n b + \) 



with the normalization constant Nf(r) being the same as that given in Eq. 161 and is plotted for differ 



ent values of 7 in Fig. 1 14 1 . 



The entropy is given by S b = - £„ fc P Ub log[P„ J and is plotted in Fig. 115. 

From Fig. fi"4) , we see that the outgoing modes have a reciprocal behaviour to that of the atomic 
system, which is expected . The superradiant behaviour of the outgoing modes is evident when one 
compares Fig. fi"4) with Fig. |[TJ, and one sees clearly an enhancement in the outgoing radiation, which 
is due to the superradiant transition of the excited atomic system formed by the quantised black hole 
modes and the incoming modes. 
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FIG. 12: Particle number distribution of the atomic modes for different times characterized by f. 
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FIG. 13: Entropy of the atomic system for different number of atomic modes characterized by j. 



V. CONCLUSION 

Laboratory models for black holes and Hawking radiation, serve two purposes. Firstly, they enable us 
to see in the laboratory, the physics of astrophysical phenomena, which are inaccessible to us exper- 
imentally, and secondly, to test various theoretical ideas associated with quantum gravity. Hawking's 
original conception of radiation from black holes, does not take into consideration the quantum struc- 
ture of the horizon of the black hole , advocated by recent work in quantum gravity 1241 , [251, [261 . A 
complete analysis of such quantum phenomena, requires a full field theoretic study of quantum gravity 
such as loop quantum gravity 124). Toy models , such as the one we have studied, give clues to how 
Hawking radiation would be modified if black hole modes are quantized. In the spirit of Beckenstein 



18 



Particle Number Distribution 
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FIG. 15: Entropy of the outgoing modes, as a function of the time characterized by y. 

and Mukhanov |25|, who made the proposal of treating black holes as quantum atomic systems and 
considering the resultant modification in the character of thermal Hawking radiation, we have tried 
to see how this can be realized in a quantum optical system that mimics such a proposal. Firstly we 
have considered the fact that a two-mode squeezed sytem of a parametric amplifier, can give rise to a 
"Hawking-like" distribution for the outgoing (signal) modes, this corresponds to the standard Hawking 
picture, where the black hole is considered to be classical. Then, we have included the quantization of 
the black hole modes to show that, indeed, the thermal character of the Hawking radiation, is modified. 
Moreover we have shown how the entanglement properties change without and with quantization of 
the black hole modes. Indeed, the figures for the entropy of the black hole modes (Figj5|and that for the 
outgoing modes shown in (Fig|9j, are suggestive of a final stable state for the black hole entangled with 
the radiation measured at infinity. 
The Hamiltonian that we have considered can be viewed in two different ways, one of which admits 
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a quadratic polynomially deformed su(l,l), while the other admits a quadratic polynomially deformed 
su(2) spectrum generating algebras. The utilization of quadratic polynomial deformations of su(l,l) 
algebras has enabled us to study the entanglement properties of the three mode quantum system (in- 
coming modes, black hole modes and the outgoing modes) in the long time limit, in an exact fashion. 
We find that the system has superradiant properties that are akin to atomic "Dicke superradiance" by 
means of a mapping of the system to a quadratic polynomially deformed su(2) algebra. 

Since, we can control quantum optical systems and atom-radiation systems (cavity systems) are easily 
created, modification effects on the Hawking radiation due to black hole quantization can be actually 
realised in this laboratory analogue. We hope that this study of a simple analogue would help in giving 
some direction to the complete field theoretic studies in quantum gravity and also help in establishing 
the similarity in the behaviour of superradiance from atomic systems and superradiance from black 
holes. In fact, in any canonical quantum gravity study, ultimately a mode decomposition in terms of 
creation and annihilation operators of the field is done, so the symmetries, in the form of deformed 
polynomial algebras, could prove to be useful for more comprehensive studies. Our ideas could also be 
applied to other type of black hole analogue systems such as those using Bose Einstein Condensates. 
Such studies, we feel, are worth pursuing, and we intend to follow these in our future work. 
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